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From  stirring  to  mixing  of  momentum:  cascades  from 
balanced  flows  to  dissipation  in  the  oceanic  interior 

James  C.  McWilliams1,  M.  Jeroen  Molemaker1,  and  Irad  Yavneh2 

Abstract.  Under  the  influences  of  stable  density  stratification  and  Earth’s 
rotation,  large-scale  flows  in  the  ocean  and  atmosphere  have  a mainly 
balanced  dynamics — sometimes  called  the  slow  manifold — in  the  sense  that 
there  are  diagnostic  hydrostatic  and  gradient-wind  balances  that  constrain 
the  fluid  acceleration.  The  nonlinear  balance  equations  are  a successful 
approximate  model  for  this  regime,  and  we  have  identified  mathematically 
explicit  limits  of  their  time  integrability.  We  hypothesize  that  these  limits 
are  indicative,  at  least  approximately,  of  the  transition  from  the  larger-scale 
regime  of  inverse  energy  cascades  of  anisotropic  flows  to  the  smaller-scale 
regimes  of  forward  energy  cascade  to  dissipation  of  more  nearly  isotropic 
flows  and  intermittently  breaking  inertia- gravity  waves.  In  the  oceans  these 
regime  transitions  occur  mostly  in  the  scale  range  of  0.1-10  km— -in  between 
the  mesoscale  and  fine-structure — where  Rossby  (Ro),  Froude  (Fr),  and 
Richardson  (Ri)  numbers  are  typically  neither  small  nor  large.  In  pursuit  of 
testing  this  hypothesis  we  have  revisited  several  classical  problems,  including 
gravitational,  centrifugal/syxnmetric,  elliptical,  barotropic,  and  baroclinic 
instabilities.  In  all  cases  we  find  definite  evidence,  albeit  still  incompletely 
understood,  of  fluid-dynamical  transitions  in  the  neighborhood  of  loss  of 
balanced  integrability. 

Introduction 

The  general  circulation  of  the  ocean  is  forced  by  sur- 
face fluxes  of  heat,  water,  and  momentum  primarily  at 
large  space  and  long  time  scales.  The  circulation  has 
comparably  large  and  long  scales,  as  well  as  important 
smaller  ones  associated  with  equatorial  zonal  and  lat- 
eral boundary  currents  and  with  the  dominant  insta- 
bility modes  at  mesoscales.  All  of  these  circulation  el- 
ements approximately  satisfy  geostrophic,  hydrostatic, 
and  incompressible  dynamical  balances. 

How  does  the  energy  dissipation  occur  for  the  gen- 
eral circulation  in  an  equilibrium  balance  with  the  gen- 
eration by  surface  fluxes?  Some  of  the  dissipation  un- 
doubtedly occurs  within  turbulent  boundary  layers  near 
the  surface  and  bottom.  Some  dissipation  also  occurs 
through  creation  of  internal  gravity  waves  by  flow  over 
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topography,  with  subsequent  wave  propagation  into  the 
interior  and  a wave-dynamical  cascade  (sometimes  in- 
volving breaking)  down  to  dissipation  at  small  scales. 
Each  of  these  routes  to  dissipation  involves  an  extrac- 
tion of  energy  from  the  circulation  near  the  vertical 
boundaries,  although  the  bulk  of  the  energy  resides 
in  the  vertical  interior.  A more  local  route  is  directly 
through  the  interior,  turbulent  cascade  dynamics  of  the 
circulation.  In  oceanic  general  circulation  models,  the 
local  route  to  dissipation  is  implied  by  the  use  of  eddy 
diffusivities  to  parameterize  this  cascade.  The  purpose 
of  this  article  is  to  examine  the  mechanism  for  the  local 
route  to  dissipation,  without  here  trying  to  assess  the 
relative  contributions  among  these  alternative  routes. 

Our  conceptual  view  of  the  mechanism  is  the  follow- 
ing: 

Large-  and  mesoscale  circulations  typically  satisfy 
a balanced  dynamics  (as  defined  below),  which 
have  little  interaction  with  the  inertia-gravity  wave 
field; 
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balanced  turbulent  cascades  are  very  inefficient  in 
energy  dissipation; 

there  are  explicitly  specifiable  limits  to  the  regime 
of  balanced  dynamics  that  are  violated  sometimes 
for  the  circulation; 

violation  of  these  limits  leads  to  energy  transfer 
to  unbalanced  motions; 

turbulent  cascades  are  much  more  efficient  in  their 
dissipation. 

In  this  view  the  important  bottleneck  in  the  local 
route  to  dissipation  is  loss  of  balance  and  its  evolution- 
ary consequences. 

Confirming  or  refuting  this  hypothesis  is  quite  chal- 
lenging, since  it  involves  the  connectedness  of  turbulent 
cascades  spanning  several  dynamical  regimes.  Even  di- 
agnosis of  the  degree  of  balance  can  be  subtle.  In  lieu  of 
making  a more  general  test  of  the  hypothesis  as  yet,  here 
we  focus  on  the  special  situation  of  the  linear  instability 
of  balanced  steady  currents  in  a rotating,  stratified  fluid 
in  relation  to  the  conditions  for  loss  of  balance.  Since 
fluid  instabilities  have  been  the  subject  of  much  prior 
research,  we  will  tell  the  story  from  both  historical  per- 
spectives of  early  discovery  and  personal  perspectives 
of  the  implications  for  the  hypothesis  above — skipping 
over  most  of  the  literature  in  between. 

Balanced  Dynamics 

The  essential  basis  for  the  approximations  of  bal- 
anced dynamics  is  velocity  anisotropy.  In  a rotating, 
stratified  fluid  with  Coriolis  frequency  / and  Brunt- 
Vaisala  frequency  N and  away  from  boundaries,  the 
evolution  from  general  initial  conditions  or  forcing  by 
the  process  called  geostrophic  (or  balanced)  adjustment 
leads  to  a local  anisotropy  with  u,v  w,  while  radi- 
ating away  inertia-gravity  waves.  Here  z and  w are  the 
coordinate  and  velocity  components  in  the  vertical  di- 
rection, antiparallel  to  gravity,  and  (x,  y)  and  (u,  v ) are 
their  horizontal  counterparts.  The  condition  for  this  to 
occur  are  that  the  Rossby  and  Froude  numbers, 

Ro  = V/fL  and  Ft  = V/NH,  (1) 

are  not  too  large  (where  V,  H,  and  L are  characteristic 
values  for  (u,v),  z,  and  (x,y)).  Under  these  conditions 
the  vertical  momentum  balance  is  approximately  hydro- 
static, 

<t>z  « b,  (2) 

(cj)  = p/p0  is  the  geopotential  function,  b = g[l-p/p0] 
is  the  buoyancy,  and  p0  is  the  mean  density);  the  diver- 
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gence  of  the  horizontal  momentum  balance  is  approxi- 
mately 

VV  « V-/VV'  + 2[^yy  -tly),  (3) 

which  is  called  gradient-wind  balance  (V  is  the  hori- 
zontal gradient  operator);  and  the  horizontal  velocity  is 
weakly  divergent  and  thus  can  be  approximately  repre- 
sented by  a streamfunction, 

u ~ —xpy  and  v « ipx.  (4) 

The  maximal  truncation  of  the  incompressible  (Boussi- 
nesq)  equations  consistent  with  these  approximations 
and  conservation  of  either  energy  (in  Cartesian  co- 
ordinates, ( x,y,z ))  or  potential  enstrophy  (in  isopyc- 
nal  coordinates,  ( X,Y,b ))  is  called  the  balance  equa- 
tions ( Lorenz  [I960];  Gent  and  McWilliams  [1984]). 
(Many  alternative  models  have  been  proposed  for  bal- 
anced dynamics;  among  the  better  ones,  their  similari- 
ties seem  more  important  than  their  differences.)  The 
balance  equations  contain  no  inertia-gravity  wave  solu- 
tions; so  they  are  often  taken  as  a dynamical-systems 
model  for  the  (advective)  slow  manifold.  They  have 
second-order  asymptotic  accuracy  as  Ro  ~ Fr  ->  0, 
whereas  the  traditional  geostrophic  balance  and  quasi- 
geostrophic  equations  have  only  first-order  accuracy 
( Gent  and  McWilliams  [1983]).  They  have  been  shown 
in  many  analyses  to  accurately  represent  the  observed 
state  and  evolution  of  large-scale  flows  in  the  atmo- 
sphere and  ocean;  for  example,  they  often  are  used  for 
initialization  of  weather  forecasts,  even  ones  for  hur- 
ricanes. An  important  aspect  of  this  accuracy  is  the 
weakness  of  inertia-gravity  wave  generation  by  balanced 
motions  when  Ro  and  Fr  are  not  large. 

The  advective  dynamics  of  balance  or  quasigeostro- 
phic  equations — called  geostrophic  turbulence — yields 
an  inverse  turbulent  cascade  of  energy  towards  larger 
scales  in  (x,  y,  z),  hence  away  from  dissipation  by  molec- 
ular viscosity  at  small  scales,  and  a forward  cascade  of 
potential  enstrophy  (i.e.,  variance  of  potential  vortic- 
ity)  to  its  dissipation  at  small  scales  ( Chamey  [1971]; 
McWilliams  et  al.  [1994]).  (This  behavior  is  analogous 
to  the  turbulence  in  a two-dimensional  fluid.)  In  the 
enstrophy  inertial  range,  Ro  and  Fr  are  invariant  as 
the  scale  decreases,  at  least  in  the  limit  of  Ro,  Fr  — > 0. 
There  is  as  yet  much  less  experience  with  balanced  tur- 
bulence at  finite  values  of  Ro  and  Fr,  but  available 
solutions  indicate  that  it  is  only  modestly  more  dissipa- 
tive of  energy  ( Yavneh  et  al.  [1997]).  It  remains  an  open 
question  how  consistently  Ro  and  Fr  avoid  increasing  in 
the  forward  cascade,  either  in  the  balance  equations  or 
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more  fundamental  fluid  dynamics:  at  small  scales  does 
the  cascade  in  balanced  turbulence  generate  inconsis- 
tencies with  its  justifying  assumptions  and  how  leaky 
is  the  slow  manifold  to  unbalanced  motions?  Never- 
theless, our  present  understanding  is  that  the  balance 
equations  do  not  provide  an  efficient  route  to  energy  dis- 
sipation away  from  boundaries:  they  imply  more  stir- 
ring than  mixing  for  momentum. 

Loss  of  Balance 

An  analysis  for  the  solvability  of  the  balance  equa- 
tions is  made  in  Yavneh  et  al.  [1997]  and  McWilliams 
et  al.  [1998].  To  be  able  to  integrate  forward  in  time 
from  a balanced  state,  several  necessary  conditions  must 
be  satisfied  everywhere  within  the  domain.  Where  these 
are  violated,  there  is  a change  of  type  of  the  partial 
differential  system  and  the  initial-  and  boundary-value 
problem  becomes  ill-posed.  For  the  balance  equations 
in  isopycnal  coordinates  (i.e.,  Gent  and  McWilliams 
[1984])  and  / > 0,  the  conditions  for  loss  of  balance 
are  the  following: 

1.  Change  of  sign  of  vertical  stratification,  N2  = 
db. 

dz ’ 

2.  Change  of  sign  of  absolute  vorticity,  A = f + 
C(z)  = f + vx  — uy  (where  the  horizontal  deriva- 
tives denoted  by  capital  letters  are  in  isentropic 
coordinates); 

3.  Change  of  sign  of  A — |S|  (where  S2  = (ux  - 
vy)2  + (vx  + uy)2  is  the  variance  of  the  strain 
rate). 

None  of  these  conditions  occurs  in  the  quasigeostrophic 
limit,  since  A,A-\S\-*f  + 0(Ro)  and  N ->  N0  + 
0{Ro),  where  N0(z)  is  the  resting-state  stratification. 
Note  the  greater  susceptibility  of  anticyclonic  regions 
(i.e.,  with  C(z)//  < 0)  in  the  second  and  third  con- 
ditions; furthermore,  note  the  greater  susceptibility  to 
the  third  condition,  since  A — |S|  < A.  The  first  and 
second  conditions  also  are  related  to  the  potential  vor- 
ticity, q = AN2.  Since  potential  vorticity  is  conserved 
on  parcels,  except  for  mixing  effects,  there  is  thus  an 
evolutionary  inhibition  for  an  unforced  flow  to  spon- 
taneously develop  a violation  of  the  first  and  second 
conditions.  However,  there  is  no  such  constraint  with 
respect  to  the  third  condition,  which  indicates  another 
sense  in  which  there  may  be  a greater  susceptibility  to 
the  third  condition. 


Instability  and  Cascade 

We  now  ask  what  happens  when  there  is  a loss  of 
balance  as  defined  above.  Obviously  any  further  in- 
tegration of  the  balance  equations  is  precluded.  So 
the  question  must  be  answered  in  a more  fundamental 
fluid  dynamics,  such  as  the  incompressible  Boussinesq 
equations,  which  have  no  restriction  on  Ro  and  Fr  for 
their  validity.  In  general,  our  expectation  is  that  some 
inertia-gravity  waves  and/or  more  nearly  isotropic  tur- 
bulence will  be  generated — instigating  some  degree  of 
enhanced  dissipation — where  there  is  a loss  of  balance. 
However,  the  efficiency  of  their  generation  is  uncertain, 
as  is  whether  the  subsequent  evolution  systematically 
departs  from  balanced  dynamics  or  relaxes  back  to- 
wards it  (e.g.,  as  a geostrophic  adjustment  or  selective 
cascade  and  dissipation  processes).  In  the  rest  of  this 
article,  we  address  the  issue  of  generation  and  subse- 
quent evolution  only  in  a very  particular  context,  viz., 
the  linear  instabilities  of  rotating,  stratified  flows  which 
are  steady,  inviscid,  balanced  solutions  of  the  Boussi- 
nesq equations.  While  this  is  far  from  the  general  cir- 
cumstances of  fluid  evolution,  it  does  provide  a cleanly 
posed  question  that  also  is  one  that  can  be  answered 
in  part  by  reference  to  the  extensive  literature  on  fluid 
instabilities. 

Instability  Types 

Now  we  attempt  to  interpret  the  known  instabilities 
in  relation  to  the  conditions  for  loss  of  balance  for  ro- 
tating, stratified  flows  where  Ro  and  Fr  are  not  large. 
It  is  probably  improvable  that  any  such  taxonomy  of 
instability  types  can  be  complete  and  unique;  however, 
after  all  the  research  that  has  gone  into  this  topic,  the 
landscape  has  become  fairly  well  mapped. 

Quasigeostrophic  Inflectional  Instability 

Consider  the  instability  of  a geostrophic  parallel 
flow  U(y,z)  with  background  stratification  Na{z)  and 
Coriolis  frequency  f(y)  in  the  quasigeostrophic  limit, 
Ro,  Fr  — ¥ 0.  Following  Rayleigh  [1880]  and  Drazin  and 
Howard  [1966]  and  ignoring  vertical-boundary  effects, 
one  can  derive  a “Rayleigh  theorem”  from  the  potential 
vorticity  equation  that  a necessary  condition  for  inviscid 
instability  of  a non-symmetric  (i.e.,  dx  ^ 0),  normal- 
mode fluctuation  is  that  the  mean  potential  vorticity 
gradient, 
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change  sign  within  the  domain.  This  is  called  an 
inflection-point  instability  since  it  involves  the  hori- 
zontal and/or  vertical  curvature  of  U ; depending  upon 
which  curvature  is  dominant,  the  instability  is  labeled 
barotropic  or  baroclinic  instability.  The  unstable  mode 
is  itself  geostrophically  balanced.  This  is  the  only  type 
of  instability  that  occurs  in  the  quasigeostrophic  limit. 
It  has  an  analytic  continuation  to  finite  Ro  and  Fr, 
where  it  can  be  expected  to  remain  balanced  over  some 
range  of  these  parameters.  Thus,  its  onset  conditions 
have  nothing  to  do  with  the  limits  of  balance,  and  it 
represents  a mode  of  stirring  within  balanced  dynam- 
ics. 

Gravitational  Instability 

The  condition  for  the  onset  of  gravitational  instabil- 
ity in  the  limit  of  vanishing  viscosity  is  N2  < 0 ( Rayleigh 
[1916];  Chandrasekhar  [1961]).  This  coincides  with  the 
first  condition  for  loss  of  balance,  and  the  mode  of  insta- 
bility is  unbalanced  (e.g.,  the  vertical  momentum  bal- 
ance is  non-hydrostatic). 

Symmetric  Centrifugal  Instability 

For  a balanced,  parallel  flow  U(y,z),  the  neces- 
sary and  sufficient  condition  for  inviscid  instability  of 
a parallel-symmetric  perturbation  is  a change  of  sign 
of  potential  vorticity  q(y,z)  ( Hoskins  [1974]).  This  co- 
incides with  the  second  condition  for  loss  of  balance. 
(This  is  a 2D  problem  rather  than  a 3D  one,  and 
McWilliams  et  al.  [1998]  show  that  the  third  condition 
is  not  germane  in  this  situation.)  For  a balanced,  ax- 
isymmetric,  azimuthal  flow  U (r,  z),  the  conditions  for 
the  inviscid  instability  of  an  axisymmetric  perturbation 
are  the  change  of  sign  of  either  the  absolute  vorticity 
A or  the  absolute  circulation  C = |/r  + U ( Rayleigh 
[1916];  Ooyama  [1966]).  McWilliams  et  al.  [1998]  show 
that  these  coincide  with  the  second  or  third  conditions, 
respectively,  for  loss  of  balance  in  this  case.  Thus,  the 
boundaries  for  onset  of  symmetric  centrifugal  instabil- 
ity, which  has  unbalanced  growing  modes,  occur  exactly 
at  the  limits  of  balanced  evolution. 

Elliptical  Instability 

The  inviscid  instability  of  the  balanced,  elliptical, 
two-dimensional,  barotropic  flow  in  an  unbounded  do- 
main, 

^(z,2/)  = l ipix2  + (3y2) , 1 > - > 0,  (6) 

z a 

was  originally  analyzed  for  / = N = 0 (Pierrehumbert 
[1986];  Bayly  [1986];  Craik  and  Criminate  [1986])  and 


(A-  ISI)/f 

Figure  1.  Growth  rate,  a//,  for  the  elliptical  flow  (6), 
maximized  over  vertical  wavenumber.  Curves  are  shown  for 
P/a  = 0.25  (dashed),  0.11  (dotted),  and  0.026  (solid).  The 
corresponding  abscissa  values  for  the  second  condition  for 
loss  of  balance  are  -0.6,  -0.8,  and  -0.95  (see  text). 


(A-ISI)/f 


Figure  2.  As  figure  1,  but  with  logarithmic  ordinate.  The 
evident  noise  is  because  of  intermittent  underestimates  of  a 
due  to  incomplete  optimization  searches  over  the  very  nar- 
row unstable  bands  in  vertical  wavenumber. 
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later  extended  to  Ro,Fr  < oo  by  Miyazaki  [1993]. 
The  unstable  modes  occur  in  bands  of  the  vertical 
wavenumber.  They  have  temporally  oscillatory  hori- 
zontal wavenumbers  and  exhibit  exponential  growth  av- 
eraged over  a wavenumber  oscillation  period. 

The  problem  was  revisited  in  McWilliams  and  Yavneh 
[1998]  from  the  perspective  of  loss  of  balance:  ellipti- 
cal instability  disappears  in  the  quasigeostrophic  limit 
and  its  onset  nearly,  but  not  precisely,  coincides  with 
the  third  condition  for  loss  of  balance  (see  Figures  1- 
2).  The  unstable  mode  is  unbalanced.  For  this  basic 
flow  the  vorticity  and  strain  rate  are  spatially  uniform 
and  unequal  in  magnitude.  The  abscissa,  (A  — |S|)//, 
= +1  in  the  quasigeostrophic  limit,  = 0 at  the  third 
condition  for  loss  of  balance,  = —(a  — (3)(a  + f3)  at 
the  second  condition  for  loss  of  balance,  and  > +1  for 
cyclonic  flows.  Cyclonic  elliptical  flows  are  stable  for 
a,  13  = 0(f). 

Taylor-Couette  and  Barotropic  Instabilities 

Consider  Taylor-Couette  flow  in  the  gap  between  two 
axisymmetric,  rotating  cylinders, 

U(r)  = Ar  + ~,  (7) 

r 

which  is  a viscous  steady  solution  commonly  studied 
in  laboratory  experiments.  The  classical  instability  for 
this  flow  is  centrifugal  ( Taylor  [1923];  Chandrasekhar 
[1961]),  whose  onset  in  the  inviscid  limit  coincides  with 
the  second  condition  for  loss  of  balance.  In  a neu- 
trally stratified  fluid,  this  is  the  only  type  of  linear  in- 
stability, since  this  profile  does  not  have  an  inflection 
point.  However,  for  a stably  stratified  fluid  with  small 
Fr,  this  barotropic  shear  flow  has  another  class  of  un- 
balanced instabilities  at  finite  Ro  in  the  anticyclonic 
regime  ( Molemaker  et  al.  [2001];  Yavneh  et  al.  [2001]), 
but  there  is  not  any  instability  for  the  quasigeostrophic 
limit,  Ro  -»  0,  nor  for  cyclonic  flows  with  Ro  = 0(1). 
For  this  new  class  an  infinite  but  countable  set  of  un- 
stable modes  exist,  which  differ  in  their  cross-stream 
structure,  each  with  a different  narrow  band  of  verti- 
cal wavenumbers.  In  Figures  3-4,  the  growth  rates  are 
shown  for  the  first  three  modes  in  the  thin-gap  limit 
where  and  S are  nearly  uniform  and  equal  in  mag- 
nitude. As  with  elliptical  flow,  the  strength  of  the  in- 
stability strongly  increases  in  the  neighborhood  of  the 
third  condition  for  loss  of  balance.  (And,  in  this  case, 
(A-|5|)//  « —1  corresponds  to  the  second  condition.) 
The  instability  can  be  shown  to  involve  a resonance  of 
shear-modified  neutral  modes  (as  required  in  any  eigen- 
value problem  whose  eigenfrequency  is  either  real  or  a 


(A — IS  l)/f 

Figure  3.  Growth  rate,  a/f,  for  Taylor-Couette  flow  (7)  in 
a narrow  gap,  maximized  over  vertical  wavenumber.  Curves 
are  shown  for  the  three  gravest  modes  in  cross-stream  struc- 
ture. 

complex-conjugate  pair);  for  the  gravest  unstable  mode 
the  resonance  involves  a pair  of  Kelvin  waves  propagat- 
ing cyclonically  along  each  cylinder  wall,  and  for  the 
other  modes,  one  of  the  Kelvin  modes  is  replaced  by 
an  inertia-gravity  mode.  Thus,  the  instability  has  an 
unbalanced  dynamics.  An  explicit  formula  can  be  ob- 
tained for  the  unstable  growth  rate  of  the  gravest  mode, 

- ~ e-i'Ro  ~ e-2^1'*),  (8) 

asymptotically  as  Ro  = max  \C<'Z')\/ f — > 0+  and 

X — (A  — |S|)//  —>  1“  (with  7 = 2 (analytically) 
for  the  gravest  mode  and  7 « 3 (computationally)  for 
the  higher  modes;  n.b.,  X is  the  abscissa  in  Figures  3- 
4).  This  asymptotic  regime  is  accurately  realized  even 
in  the  neighborhood  of  the  third  condition  for  loss  of 
balance,  which  implies  that  there  is  an  extremely  rapid 
weakening  of  o with  Ro  and  X in  this  neighborhood, 
but  not  an  abrupt  cessation  at  any  critical  value  near 
X = 0. 

Recently  we  have  also  solved  the  linear,  inviscid  sta- 
bility problem  for  a barotropic  boundary  current  along 
a straight  coastline  (at  x = 0), 

V(x ) = V0exp(-ax),  (9) 

in  a uniformly  rotating,  stratified  fluid  in  a semi-infinite 
domain  (i.e.,  x > 0).  This  profile  also  does  not  have  an 
inflection  point.  In  addition  to  no  normal  flow  at  the 
boundary,  we  prescribe  a radiation  condition  at  a suffi- 
ciently distant  location  in  the  interior  by  matching  the 
solution  to  an  outwardly  radiating,  free  inertia-gravity 


64 


McWilliams,  molemaker,  and  yavneh 


(A  — fS  I )/f 


Figure  4.  As  figure  3,  but  with  logarithmic  ordinate 
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Figure  5.  Growth  rate,  cr/f,  for  the  parallel  boundary 
current  (9),  maximized  over  along-stream  and  vertical 
wavenumbers.  The  abscissa  value  is  based  upon  the 
minimum  over  x which  occurs  at  x = 0. 


(A-  IS I)/ f 


Figure  6.  As  figure  5,  but  with  logarithmic  ordinate 


wave, 

=  •  **“>  (10) 

where  k is  a cross-stream  wavenumber  determined  from 
vertical  and  along-stream  wavenumbers  and  eigenfre- 
quency  using  the  dispersion  relation  of  an  inertia-gravity 
wave.  Again  there  are  imbalanced  unstable  modes  for 
anticyclonic  flows  (i.e.,  V0/f  >0)  away  from  the  quasi- 
geostrophic  limit.  In  contrast  with  Taylor-Couette  flow, 
with  its  discrete  spectrum  of  modes  in  the  cross-flow 
direction,  here  there  is  a continuous  set  of  unstable 
modes  for  any  along-stream  and  vertical  wavenumber 
pair  (and  its  corresponding  cross-stream  k).  In  Fig- 
ures 5-6,  growth  rates  are  shown  for  optimal  along-flow 
and  vertical  wavenumbers,  indicating  yet  again  a rapid 
increase  of  the  growth  rate  in  the  vicinity  of  the  third 
condition  for  loss  of  balance;  the  functional  form  of 
cr\Ro]  here  appears  to  be  close  to  that  for  the  higher- 
mode, Taylor-Couette  instabilities.  Again,  a resonance 
can  be  diagnosed  with  a shear-modified  Kelvin  mode 
and  an  inertia-gravity  mode. 

In  summary,  the  three  different  barotropic,  anticy- 
clonic, rotating,  stratified,  shear  flows  analyzed  in  this 
and  the  preceding  section  all  have  an  unbalanced  insta- 
bility whose  strength  rapidly  increases  in  the  vicinity 
of  the  third  condition  for  loss  of  balance  (in  addition 
to  the  imbalanced,  centrifugal  instability  in  the  vicin- 
ity of  the  second  condition).  In  the  hydrostatic  limit, 
this  instability  is  related  to  ageostrophic,  parallel,  anti- 
cyclonic flow  instability  in  a shallow-water  layer  ( Sato - 
mura  [1981];  Griffiths  et  al.  [1982],  et  seq.) 

Baroclinic  and  Kelvin-Helmholtz  Instabilities 

Consider  a steady,  spatially  uniform  vertical  shear 
flow,  U(z)  oc  z,  in  a uniformly  stratified  and  rotating 
fluid  in  a vertically  bounded  domain  with  non-isopycnal 
boundaries.  We  define  the  Richardson  number  for  this 
flow  by  Ri  = N2/U2,  which  is  equivalent  to  a Fr~2. 
The  balancing  buoyancy  field  is  B(y,  z)  = N2z  - fUzy, 
and  the  associated  quantities  in  the  second  and  third 
conditions  for  loss  of  balance  are 

*-4 -s]'  x-s  = 4“jI]'  (u) 

Notice  that  this  flow  regime  is  anticyclonic  since  A/ f < 
1.  Thus,  we  can  make  the  following  categorization  for 
this  flow  in  terms  of  Ri  and  in  relation  to  the  conditions 
for  loss  of  balance: 

• The  quasigeostrophic  limit  occurs  for  Ri  — >■  oo.  It 

has  a geostrophic  baroclinic  instability  ( Chamey 


FROM  BALANCED  FLOWS  TO  INTERIOR  DISSIPATION 


65 


[1947];  Eady  [1949]),  whose  “inflection  point”  in 
this  particular  case  occurs  at  the  vertical  bound- 
aries (see  above). 

• The  third  condition  is  satisfied  if  Ri  < 2. 

• The  second  condition  is  satisfied  if  Ri  < 1. 

• The  classical  (non-rotating)  condition  for  the  on- 
set of  Kelvin-Helmholtz  instability  ( Miles  [1961]; 
Howard  [1961])  is  satisfied  if  Ri  < 

• The  first  condition  is  satisfied  if  Ri  < 0;  this  is 
gravitational  instability  (see  above). 

Analyses  of  this  problem  in  Stone,  [1966,1970]  and 
Nakamura  [1988]  show  that  a centrifugal  instability  does 
occur  near  Ri  = 1 for  both  zonally  symmetric  and 
asymmetric  fluctuations  and  that  another  ageostrophic 
instability  (with  shorter  zonal  wavelengths  than  the 
geostrophic  instability)  occurs  for  even  larger  values  of 
Ri,  though  its  onset  value  is  not  well  determined.  This 
latter  instability  is  shown  in  Figure  2 of  Stone  [1970], 
with  an  accompanying  remark  that  it  has  “growth  rates 
[Im<r]  which  decrease  as  Ri  increases,  and  may  in  fact 
still  be  unstable  for  Ri  > 2,  but  if  so  the  growth  rates 
were  too  small  to  to  be  found  by  our  numerical  method, 
because  of  the  near  singular  behavior  of  the  coefficients 
of  the  differential  equation  when  Imcr  is  very  small.” 
Confirmation  is  presented  in  Figure  8 of  Nakamura 
[1988]  which  shows  that  Imu  is  strongly  decreasing  as 
Ri  increases;  it  is  somewhat  unclear  exactly  how  large 
a Ri  value  he  obtained  unstable  solutions  for,  but  his 
figure  suggests  it  is  at  least  as  large  as  Ri  = 2.  An 
interpretation  of  the  shortwave  instability  is  presented 
(p.  3261):  it  involves  a resonance  between  a boundary- 
trapped  shear  mode  and  an  inertia-gravity  mode,  with 
an  “inertia  critical  level”  that  limits  the  vertical  extent 
of  the  eigenmode.  Thus,  although  further  examination 
is  needed  of  the  behavior  of  Imcr  in  the  vicinity  of  the 
third  condition  for  loss  of  balance,  the  instabilities  for 
this  flow  appear  to  support  quite  well  the  hypothesis  we 
have  advanced. 

The  Rest  of  the  Quest 

We  have  presented  here  a survey  of  the  known  insta- 
bilities of  rotating,  stably  stratified  fluids.  Barotropic 
and  baroclinic  inflectional  instabilities  are  the  only 
types  that  occur  in  the  quasigeostrophic  limit,  and 
several  others  arise  in  the  neighborhood  of  the  limits 
of  balance  as  determined  from  the  balance  equations. 
These  latter  include  gravitational,  symmetric  centrifu- 
gal, and  elliptical,  as  well  as  presently  unnamed  variants 


of  barotropic  and  baroclinic  instabilities.  The  insta- 
bility onset  is  generally  understood  to  be  sharp  with 
respect  to  the  N2  and  A conditions  for  loss  of  bal- 
ance, which  correspond  to  the  gravitational  and  cen- 
trifugal instabilities.  We  have  shown  for  the  elliptical, 
Taylor-Couette,  and  barotropic  boundary-current  insta- 
bilities that  the  onset  behaviors  are  smooth  in  the  tran- 
sition across  the  A — (S'!  condition,  though  still  quite 
steep,  with  exponential  or  steeper  dependences  for  the 
growth  rate,  a(A  — |S|).  The  onset  behavior  near  the 
A—  |S|  condition  has  yet  to  be  as  well  determined  in  the 
ageostrophic  baroclinic  instability  problem,  but  previ- 
ous studies  indicate  it  is  also  steep.  Kelvin-Helmholtz 
instability  does  not  fit  in  this  categorization  scheme  be- 
cause its  onset  occurs  well  beyond  the  A and  A — |S| 
conditions. 

Thus,  we  conclude  that  these  problems  give  confir- 
mation, pro  tern,  of  the  hypothesis  that  the  limits  of 
balance  are  indicative  of  transition  from  the  larger-scale 
regime  of  inverse  energy  cascades  in  anisotropic  flows  to 
the  smaller-scale  regimes  of  forward  energy  cascade  to 
dissipation  of  more  nearly  isotropic  flows  and  intermit- 
tently breaking  inertia-gravity  waves.  There  is  still  need 
for  refining  our  understanding  of  a{Ri,  Ro,  Ft ) in  the 
neighborhood  of  the  conditions  for  loss  of  balance,  espe- 
cially for  the  less  familiar  third  condition.  Now  that  we 
know  that  pure  barotropic  and  baroclinic  instabilities  fit 
the  hypothesis,  it  would  be  worthwhile  to  further  test 
it  with  more  general  shear  profiles  (e.g.,  the  “coastal” 
V (x,  z ) profile  that  Barth  [1994]  has  shown  to  have  an 
ageostrophic  baroclinic  instability).  Beyond  these  lin- 
ear instability  problems  lies  the  challenge  of  aptly  diag- 
nosing flow  evolution  near  the  loss  of  balance  conditions 
in  general  nonlinear  initial-  and  boundary-value  prob- 
lems. Nevertheless,  we  interpret  the  results  in  hand  as 
indicating  that  the  loss-of-balance  transitions  are  steep 
but  fuzzy,  consistent  with  the  view  that  the  balanced 
slow  manifold  is  itself  modestly  fuzzy  and  leaky  to  more 
efficiently  cascading,  unbalanced  motions. 
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